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FOREWORD 

The  research  described  in  this  report  was  performed  in  fiscal  year 
1986  as  a  part  of  an  effort  to  apply  robust  statistical  estimation 
techniques  tc  radar  inverse  scattering.  The  goal  in  radar  inverse 
scattering  is  to  identify  an  unknown  radar  target  by  techniques  that 
make  fundamental  use  of  the  physics  of  electromagnetic  scattering. 
Inverse  scattering  raethoas  yield  basic  descriptor  of  targets,  such  as 
lengths,  diameters,  and  radii  of  curvature.  The  hope  is  that  these 
descriptors  will  permit  target  identification  without  the  computer 
intensive  signal  template  matching  required  by  conventional  target 
identification  methods.  This  effort  was  supported  by  6.1  funds  from 
Naval  Air  systems  Command,  Code  933;  support  for  this  work  has  con¬ 
tinued  in  fiscal  year  1987  under  funding  fror  the  Office  of  Naval 
Research,  Code  1247. 
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INTRODUCTION 


Many  approaches  to  electromagnetic 
reconstruction  of  target  shape  only  by  assuming  broadband  data  over  a 
large,  or  even  complect-,  set  of  aspects  (References  1  and  2).  Because 
of  problems  with  non uniqueness  of  solutions,  this  is  often  necessary. 
Unfortunately,  howev.tr,  such  stringent  requirements  limit  the  applica¬ 
bility  of  the  developed  techniques.  In  radar-target-identification 
problems,  for  example,  we  typically  are  restricted  to  high  frequencies 
and  an  unknown  but  ve*y  small  range  of  aspects  (aspects  being  defined 
by  target  orientatic  ij.  Adding  to  these  difficulties  are  typical 
complex  target  structures  consisting  of  distributed  (perhaps  interac¬ 
ting)  scatterers  v'  c’.i  varying  reflectivities  and  associated  intrinsic 
phase  shifts.  In  .'net,  very  few  target  characteristics  can  be  deter¬ 
mined  with  certainly;  we  must  approach  the  problem  statistically. 


Surprisingly,  statistical  methods  in  theoretical  inverse  scatter¬ 
ing  have  not  found  much  favor  in  problems  involving  target-shape  esti¬ 
mation.  Moreover,  many  of  the  applications-oriented  developments  tend 
to  be  ad  hoc  in  nature  and  rely  heavily  on  large  data  bases  and  the 
special  training  of  human  operators.  These  methods  are  often  based  on 
the  amplitude  statistics  of  the  scattering  centers  and  frequently 
require  extensive  a_  priori  information  before  they  can  be  considered 
useful. 


We  develop  a  statistical  inverse  method  that  allows  for  shape 
reconstruction  of  such  symmetric  complex  targets.  These  targets  are 
assumed  t...  consist  of  many  independent  scattering  centers  distributed 
within  a  finite  region  of  space  and  are  expected  to  "wobble"  (with 
time)  in  a  small  but  unknown  manner.  The  data  used  in  the  reconstruc¬ 
tion  consist  of  phase  derivatives  of  the  kind  that  are  currently  used 
fer  phose-dif ference  radar  tracking.  We  proceed  by  first  examining  the 
statistics  of  these  phase  derivatives,  then  showing  how  they  are 
related  to  target  shape,  and  finally  by  developing  a  shape- 
reconstruction  algorithm  based  on  these  results.  Examples  using 
synthetic  data  are  displayed.  We  finish  by  describing  a  proposed 
method  for  measuring  the  required  phase-derivative  statistics  with  a 
practical  antenna  array. 
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STATISTICS  OF  THE  PHASE  DERIVATIVES 


A  statistical  model  based  on  cross  sections  can  hatdly  be  expected 
to  be  of  much  utility  in  shape  determination  tor  complex  targets, 
because  scattering  amplitude  depends  in  a  sensitive  and  complicated  way 
on  local  scattering  center  structure  and  morphology.  However,  it  has 
long  been  realized  that  a  wealth  of  target-shape  information  is 
contained  in  the  phase  of  the  asymptotic  scattered  field.  The  trouble 
is,  of  course,  that  this  phase  information  is  very  difficult  to  obtain 
with  any  accuracy  over  a  set  of  aspects.  On  the  other  hand,  phase- 
derivative  information  can  be  approximated  to  good  accuracy  by 
measuring  the  phase  differences  between  closely  spaced  sensors.  In 
radar  applications,  for  example,  the  method  is  known  as  phase  roonopulse 
and  hjs  been  extensively  studied  in  its  relationship  to  tracking 
problems  (Reference  3).  Because  of  practical  considerations,  such  as 
noisy  measurements  and  inaccurate  aspect  information,  the  first 
derivatives  are  determined  by  these  kinds  of  "point"  measurements  using 
two  receiving  antennas.  In  general,  the  derivative  would  require 

an  array  of  J  +  1  antennas. 

If  F  =  denotes  the  scattered  field,  then  the  phase  is  just 

♦  *  arctan  (-1  TT'f) 


so  that 


1£  -  F3F*/3n  -  F*3F/3n 

3n  “  1  2FF* 


Consider  a  situation  in  which  the  scattering  centers  associated  with  a 
complex  extended  body  are  excited  by  an  incident  plane  wave  originating 
from  an  observer  located  at  a  distance  R  from  the  body.  (The 
coordinate  system  is  fixed  within  the  scatteier. )  The  high-frequency 
time  harmonic  far-field  scattered  from  a  general  target  can  be  written 
in  the  form  (Reference  4) 

E  ei2kR  N  i(2kR  -  r.  +  4>. ) 

F<“>  ’  -tst-  *  V  J  J 

j-1  J 
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where  the  scatterer  is  taken  to  be  composed  of  N  (possibly  infinite) 
scattering  centers  with  associated  (real)  scattering  strengths  Aj , 
positions  rj,  and  intrinsic  phase  shifts  4>j.  E0  is  the 

(constant)  magnitude  of  the  incident  field,  R  is  a  unit  vector  directed 
away  from  the  observer,  k  «  2tt /X  is  the  wave  number,  and  2(R  •  ri  + 
R)  is  the  (two-way)  wave-travel  distance  from  the  observer  to  the  j^h 
scatterer.  (We  have  suppressed  the  e^ll)t  time  dependence.) 

In  a  standard  cartesian  coordinate  system  in  which  R  lies  along 
the  z  direction  and  rj  has  components  (uj ,  vj ,  w-j),  define 


E 

mi  ,  in  2  j  ®  3 


i<t> 

e  ml>n2>m3 


■  V? 


i( 2kR  - 

tio  m 
v  *  w  J  P 

J  J 


+  V 


where  the  iu  are  integers.  In  terms  of  these  definitions,  the 
various  phase  derivatives  can  be  determined  from  Equation  1.  In  the 
high-frequency  limit,  we  have  for  example: 


Xioo 


X200 


X101 


3x  \  E  / 


cos(  4»  1 00) 


sin(  d>~ d>2  0  0^ 


32<j) 
3x  3k 


/E  1 0  l\ 
\  E  / 


sin(  $>--)>  10  j) 


00) 

'  (lO(^~)slnU~'hoo) 


(2) 


1 


and  so  forth.  (We  assume  that  only  the  centra]  element  of  the  antenna 
array  is  u3ed  for  estimating  wave-number  derivatives.) 

To  obtain  the  statistics  of  the  X  £mn  *  we  ma^e  t'lfc  following 
assumptions.  We  take  (Aj }  to  be  a  set  of  independent  random 
variables,  each  Aj  having  tne  same  statistical  frequency  distribution 
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I 


v. 
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f  A^Aj ) •  For  our  complex  physical  scatterer,  we  also  take  {$j  }  to 
be  a  set  of  independent  random  variables  symmetrically  distributed  on 
( — it ,  it )  with  distribution  £ ( 4>i ) •  This  assumption  is  not  only 

appropriate,  it  greatly  simplifies  the  following  calculation. 
Moreover,  at  high  frequencies  (for  which  X  is  much  less  than  a 
"typical"  target  dimension),  it  clearly  will  not  alter  significantly 
the  results  that  would  be  obtained  by  assuming  the  target  to  be 
composed  of  noninteracting,  convex,  and  perfectly  conducting  scattering 
centers.  Finally,  we  suppose  that  the  aspect  angle  through  which  the 
target  rotates  is  sufficiently  small  so  that  the  statistics  are 
stationary  and  so  that  each  A-i(R)  and  4m(R)  remain  essentially 
fixed. 

We  begin  by  examining  the  statistics  of 


cos(4>  -  A 


Define 


=  ljAj  cos(Bj),  t  =  ^Aj  oin(8j) 


u  5  cos(  ) ,  and  v  =  sin(  ) 


where  ctj  =  (2kuj /R)  ^(2kvj  /R)m( 2wj  )n,  0j  =  2kR  -  rj  + 
4>j,  and  the  sums  are  taken  to  range  from  1  to  N.  (Suppressing  the 
(<i,m,n)  dependence  will  not  cause  confusion  in  the  following.)  Then  we 
r.an  wr  Lee 


su  +  tv 
“  s2  +  t2 


Under  the  stated  assumptions,  we  can  show 


<s>  -  <t>  -  <u>  -  <v>  -  0 


o2  =  <s2>  -  <s>2  ■  a~,  and  a2  "  a2 
s  t  u  v 


vv.-vv;.  -.  -  - 


It  is  easy  to  see  (by  calculating  cross  correlations)  that  s,  t,  u,  and 
v  are  statistically  independent.  Since  each  of  these  terms  can  be 
written  as  a  sum  S  *  E;;j,  where  c-j  are  identically  distributed  with 
nonvanishing  finite  variances,  they  each  satisfy  the  Lindberg 
conditions  (Reference  5)  and  so  have  asymptotically  normal  frequency 
distributions. 

It  is  a  straightforward  calculation  to  show  that  both  E  and 
of  Equation  3  are  independent  Rayleigh-distributed  random  variables. 
The  independence  of  $  and  follows  from  the  independence  of  the 

normally  distributed  complex  variables  E  exp(i$)  and  E£,mn 

expCicJijnjn) ,  so  that  the  product  cos($  -  normally 

distributed  and  independent  of  E.  Finally,  we  can  conclude  that  the 
statistics  of  5jjjnn>  a£!  the  quotient  of  normal  and  Rayleigh- 

distributed  independent  random  variables,  can  be  expressed  as  a 
Student's  "t"  distribution: 


where 


gton(^ 


SI 


ton 


<S1L  +  «2>3/2’ 


W 


ton 


u 

2 


To  write  this  in  a  more  tractable  form,  observe  that  a  scatterer  of 
strength  M  x  A,  located  at  r j ,  has  scattered  intensity  equivalent  to 
M  (in-phase)  scatterers  of  strength  A  at  the  same  location.  Define  a 
weighted  density  function  by 


A*p(r  )AV  =  /  A2  £A(A  )dA 


where  A  is  a  constant  and  AVj  is  an  incremental  volume  about  rj. 
Substitution  yields 
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lim 

N+“> 


A*  li  avi 

A2  IjPCjMVj 


/"  a2  p(r)d3r 


(See  Reference  5 


for  the  details  of  these 


calculations. ) 


Of  course,  this  result  gives  the  frequency  distribution  of  the 
first  derivatives  of  the  phase.  The  frequency  distributions  of  the 
higher  derivatives  follow  from  this  result.  for  example,  from 
Equation  2  we  see  that  the  second  derivative  of  with  respect  to  x  can 
be  treated  as  a  function  of  three  independent  Student's  "t"  distributed 
random  variables : 


z  -  sx  -  s2s3 


We  can  compute  the  statistical  distribution  of  z  in  the  usual  way  as 


00  oo 

f  (z)  -  f  r  fl  oo  a  (Z  +  r)  a,nn(~)  «inn(t)dtdr 

y  *  *  —  v  u  ^  *  w  u  ^  —  a  w  v 

—  OO  —  OO 


where  the  gj^jm  are  defined  by  Equation  4.  Making  the  substitution 
n  r  t2/r,  the  integral  over  t  can  be  seen  to  be  of  the  tabulated 
variety  (Reference  6).  After  some  algebra,  we  are.  left  with 


f 

z 


fl2200  B(3/2,3/2) 
2[(n2lOQ-z'2  +  0  0^ 


f1  lL_--y->- 
i-i  (1  +  y) 


g(  1/2,3/ 2;  2;  y2)  dy 
(y2  +  2by  +  a) 2 


where 
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c  n2 100  +  *)2  +  n2200 

a  =  - - - 

C  fi2 100  -  z)  2  +  fi2200 

( fl2 100  -  z)^2100  +  z'  -  ^2200 
b  =  - 

(fl2 100  ~  z)  2  +  fl2200 

r  ”  sj2ioo 

y  =  - - - 

r  +  Q  ioo 

is  a  beta  function,  and  F  is  a  hypergeometric  function.  We  shall 
satisfy  ourselves  with  an  approximation  to  this  last  Integral. 

Observing  that  the  term  (1  -  y)3F(l/2,  3/2;  2;  y 2 ) / ( 1  +  y)  is 
singular  at  y  *  -1  and  decays  very  rapidly  to  0  thereafter,  we  take 


f  «  C  j1  S(y  +  lKy2  +  2by  +  a)~3/ 2  dy 
z  —  1 


where  C  is  a  normalization  constant  and  6  is  a  delta  function.  Our 
final  result  follows  easily  from  this  last  integral. 

We  have  shown  that  the  distributions  of  the  phase  derivatives 
Xgmn  (F.quation  2)  obey 


f  100<x)  •  J  «ioo2^\oo2  +  x2>“3/2 


f200^X)  “4  n200£(R2C02  +  X2)'3/2 


fioi^x)  “7  ^ioi2^ioi2  +  x2)-3/2 


(5) 


etc.  ...,  where 
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«1  OO2  ■  (f-)  /  u2  P<r)d3r 

«2002  “  /  u4  P(r)d3r 

«1 0 1 2  “  22  /  u2w2  p(r)d  ?r 

The  distributions  of  the  higher  derivatives  follow  from  similar 
considerations. 

In  the  next  section,  we  shall  show  how  the  can  be  related 

Jtmn 

to  target  structure. 


TARGET-SHAPE  DESCRIPTION 

In  the  following,  we  take  as  cur  target— shape 
scattering  center  density  function  p(r).  Associated 
characteristic  function  defined  by 


9(  A)  -  /  p(r)e1^  r  d3r  (6) 


^  ~  —  *-U  ~ 

U1IC 

with  this  is  a 


Expanding  9(A)  by  Maclaurin's  formula  yields 

6(A)  -  1  +  l  in  l 

n« 1  m  i  ,ra2>ni 

where 

=  ±-n  f  3n9(A) 

ym1,m2,m3  =  l3A1ml9A<;ni23A3m3 


A  -  0 


,raiA,m2A/3 


+  Rv  (7) 
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0.  Subotitui.  .'.ou  ol  Equation 


and  R.,  is  a  r  ■•painder  tei  ir:  lin  R,,  • 

v  y-Mn  v 

6  into  this  last  definition  yields 


li 


iij,ra2,ni3 


j  um lvm2w'J3  p(r)d3r 
-  00 


(For  example,  the  moments  of  o(r)  determine  p(r)  by  Inversion  of  the 
transform  Equation  6.) 

Ta-gets  symmetric  across  the  coordinate  planes  will  have  vanishing 
"odd"  moments.  Sach  targets  can,  therefore,  be  reconstructed  from  a 
knowledge  of  their  ’even"  moments  only.  We  have  seen  that  the  statis¬ 
tics  of  the  ;  uo.se  derivatives  depend  in  a  simple  way  on  these  even 
moment'..,  and  it  is  a  straightforward  matter  to  extract  this  information 
from  a  phnse-derl vative  data  set. 


ESTIMATION  OF  P.2 

turn 


Given  a  collection  of  M-phase-derivafive  data  points  gathered  from 
the  target  modeled  in  the  previous  discussion,  we  seek  to  estimate 
nJlmn2  *  (2k/R)2li+m)2n  y2i,2ra,2n*  Ft>r  a  complex  target,  it 

is  known  that  x$mn  i-8  a  rapidly  varying  function  of  aspect.  There¬ 
fore,  we  are  often  in  the  position  of  looking  at  a  widely  varying  data 
set  from  a  target,  which  itself  varies  only  slightly  in  presented 
aspect.  For  this  situation,  we  can  determine  the  value  of  Sl^n  that 
best  matches  these  data  to  the  Student's  "t"  distribution  law  that  is 
expected. 

Assuming  the  statistics  of  Equation  5,  the  maximum-likelihood 
estimator  for  Q^un  will  be  the  one  that  maximizes  the  function  (Equa¬ 
tion  7) 


MSI,  ) 
Ism 


M 

R 

i-1 


n‘  /  2 

tan 


(n 


tan 


*i2> 


3/ 2 


where  the  product  is  taken  over  all  the  M  elements  of  the  data  set 
Here,  we  choose  instead  the  value  of  figjnn  that  maximizes 
the  logarithm  of  A( Qgmn )  (which  will  also  maximize  Q^^),  for 
obvious  reasons. 
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Differentiating  £nA(fl^inn)  with  respect  to  and  setting  the 

result  to  zero  yieids  the  relation 


2 


M 


l 

i=l 


ton 


,  3/ 2 


ton 


M 

■x 


Solving  this  equation  is  tantamount  to  finding  the  real  roots  of  a 
polynomial  cf  degree  2M,  but  since  it  Is  the  sum  of  functions  of  the 
form 


n2/2 
A2  +  x2 


1_ 

3 


then  there  will  be  only  two  real  zeros,  which  will  be  located  symmetri¬ 
cally  about  the  origin,  and  the  magnitude  of  either  zero  is  the  solu¬ 
tion.  Thus,  we  can  use  an  iterative  scheme  to  solve  for 


SAMPLE  RESULTS 


Targets  of  the  kind  that  we  hava  been  modeling  allow  for  the  easy 
construction  of  simulated  data  sets.  This  is  because  we  are  not  trying 
to  locate  the  actual  scattering  centers  associated  with  a  specific  tar¬ 
get.  Instead,  we  are  free  to  create  structured  collections  of  scatter- 
ers  and  investigate  tho  ability  of  the  devised  algorithm  to  accurately 


n V*  ao  f 
uuw  w 


plish  such  a  verification,  we  developed  a  computer  simularlon  that  (1) 
creates  a  statistical  target  model,  (2)  determines  the  relevant  phase- 
derivative  information  over  a  small  range  of  random  aspects,  (3)  esti¬ 
mates  the  corresponding  moments  of  the  scattering  center  distribution 
function,  and  (4)  reconstructs  the  "target"  by  inverting  Equation  6. 


Each  of  the  targets  displayed  in  Figures  1  through  3  consisted  of 
100  scattering  centers  randomly  located  within  a  plane-symmetric 
support  "shape.”  Each  scattering  center  was  randomly  assigned  a  scat¬ 
tering  strength  and  associated  phase  shift.  Data  sets  consisting  of 
the  first,  second,  a.d  third  derivatives  of  the  scattered  field  phase 
were  constructed  by  "viewing"  the  target  from  different  aspects.  We 
have  chosen  these  aspects  to  lie  within  a  small  "cone"  of  directions 
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(»  3  x  3  degrees)  centered  on  one  of  the  multiple  symmetry  directions 
of  the  target. 


(a) 


FIGURE  1.  Synthetic  Target  (a)  and  Its  Reconstructed  Contour  "Image" 
(b).  The  target  consisted  of  100  scattering  centers  randomly  placed 
inside  the  displayed  support  with  random  scattering  strengths  (shown  in 
the  figure)  and  random  local  phase  shifts.  Target  dimensions  are  10  x 
10  x  1;  modeled  wavelength  is  X  ■  10-2;  and  range  ir  R  ■  10 3.  The  data 
consisted  of  the  first  three  phase  derivatives  "collected"  randomly 
from  within  a  3-  x  3-degree  range  of  aspects. 
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FIGURE  2.  Synthetic  100-Point  Spherical  Target  (a)  of  radius  10  and 
Its  Reconstructed  Image  (b).  All  other  parameters  are  as  ’.v  Figure  1. 
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FIGURE  3.  Rectangular  Plate  Target  (a)  and  the  Image  Recon¬ 
structed  From  Its  Synthetically  Determined  Scattering  Data. 
The  support  dimensions  of  the  target  are  15  x  20  x  1;  all 
other  parameters  are  as  in  Figure  1. 
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The  target  images  were  recaptured  by  using  Equation  7  to  load  a 
(64  x  64  x  64)  array  and  then  using  a  three-dimensional  fast-Fourier 
transform  algorithm  to  invert  Equation  6.  Finally,  the  recaptured 
density  function  was  collapsed  along  an  axis  to  facilitate  the  plotting 
routine;  thus  the  figures  really  represent  an  average  density  along 
this  axis. 


MEASUREMENTS  OF  THE  PHASE-FRONT  DERIVATIVES 


For  Implementation  in  a  practical  radar,  the  technique  presented 
in  this  report  requires  the  measurement  of  the  first  three  or  so 
spatial  derivatives  of  the  echo  phase  front.  This  section  discusses  a 
technique  for  measuring  these  derivatives. 


To  measure  the  phase-front  derivatives,  the  phase  front  is  sampled 
and  the  finite  difference  approximation  to  the  derivatives  formed.  We 
assume  for  this  discussion  that  the  first  three  derivatives  are 
desired.  Measurement  of  three  derivatives  requires  that  the  phase 
front  be  sampled  by  four  antennas,  as  shown  in  Figure  4,  where  for  now 
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FIGURE  4.  Four-Element  Antenna  Array  With  Impinging  Wave. 
A  sample  locus  of  constant  phase  is  shown.  The  ^ 
represent  the  electrical  distance  to  the  conscarit 

phase  locus,  where  is  the  physical  distance. 
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to  the  direction  to  the  target.  That  is,  the  tracking-angle  error  is 
zero.  Denoting  the  phase  at  each  antenna  as  (Jn°  (i  *  1...4),  the 
finite  difference  approximations  to  the  first  three  differentia.’ a  are 
given  by 

(t3) 

A$U)  -  *4  -  2^3  +  *2  (9) 

and 

A$^^  “  ^4  “  3$  3  +  3(^2  —  ^*1  (10) 

These  differentials  can  be  measured  by  the  scheme  showr  in  Figure  5. 
The  phase  discriminators  produce  an  output  that  is  the  radio  frequency 
(RF)  phase  difference  between  the  two  adjacent  antenna  elements.  The 
path  lengths  from  antenna  elements  to  discriminators  are  constructed  to 
be  equal.  The  arrangement  of  difference  amplifiers  (whose  bandwidth 
need  only  accommodate  the  rate  of  change  of  the  RF  phase,  not  the  RF 
itself)  forms  Equations  8  through  10. 


FIGURE  5.  Interconnection  of  Four  Antenna  Elements  to  Extract  the 
First  Three  Phase-Front  Derivatives. 
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In  general,  some  amount  of  tracking  error  will  always  exist 
(Figure  6),  in  which  case  (in  the  far  field)  the  RF  phase  at  each 
antenna  is  given  by 


,  .0 

?1  ”  91 

(j»2  “  ^2  ~  W  8ia  0 

m  $3  —  2  0d  sin  0 
$4  -  “  3  0d  sin  0 


FIGURE  6.  The  Phase  Front  of  Figure  4  Incident  at  an  Angle 
0.  The  squares  present  the  new  phase-front  sampling  points; 
the  triangles  are  the  sampling  points  in  Figure  4. 
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The  zero  reference  for  the  phase  is  taken  arbitrarily  at  element  1. 
Note  that  for  an  identical  wave  front,  different  values  of  the  4>^° 
are  measured,  because  rotation  of  the  wave  front  through  an  angle  0  has 
the  effect  of  decreasing  the  spacing  between  the  antennas.  This  causes 
the  wave  front  to  be  sampled  at  a  different  set  of  points.  The 
<j)g0,  however,  still  represent  the  required  phase  values  for  forming 
the  derivatives.  Forming  Equations  8  through  10,  Equation  8  becomes 


*°3 


♦s 


-  8d  sin  0 


(ID 


and  Equations  9  and  10  are  unchanged.  Hence,  a  biased  estimate  of 
AiJ>(  1 )  results,  h?!t  £<p(^)  and  £^(3)  are  unbiased.  However,  forc¬ 
ing  Equations  8  through  10  to  have  zero  mean  (necessary  conditions  for 
this  technique)  will,  for  constant  0,  remove  the  bias  term  in  Equation 
11.  Clearly,  problems  may  arise  if  the  tracking  loop  allows  fluctua¬ 
tions  in  0  that  approach  the  rate  of  change  of  the  phase  differentials, 
thus  causing  0  to  vary  during  the  measurement  interval.  The  tracking 
loop  must  be  designed  to  prevent  this  in  an  application  of  the  tech¬ 
nique  described  in  this  report. 


CONCLUSION 


We  have  developed  a  statistical  high-frequency  inverse  scattering 
method  that  allows  the  reconstruction  of  very  complex  targets  from 

extemely  limited  data.  These  data  consist  of  various  phase  derivatives 
collected  over  a  small  range  of  aspects  and  frequencies.  Signifi¬ 
cantly,  we  do  not  require  that  this  aspect  dependence  be  known.  We 
have  assumed  tne  target  to  be  composed  of  many  individual  scattering 

centers,  each  with  potentially  different  morphological  properties; 
however,  our  result  does  not  require  that  any  specific  assumptions  be 
made  about  the  statistics  of  the  local  scatterer  strengths. 

Because  we  have  used  the  classical  central-limit  theorem  to  obtain 
the  relevant  statistics  (and  from  these  the  even  moments  of  the 

scattering  center  distribution  functions),  the  technique  allows  only 

the  reconstruction  of  symmetric  targets. 
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